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Abstract 

This paper presents a realistic simulation of 
an aircraft in flight using the AD 100 digi- 
tal computer. We discuss specifically the im- 
plementation of three model features: (1) a 
large aerodynamic data base (130,000 func- 
tion values) which is evaluated using function 
interpolation to obtain the aerodynamic co- 
efficients, (2) an option to trim the aircraft 
in longitudinal flight, and (3) a flight con- 
trol system which includes a digital controller. 
Since the model includes a digital controller 
the simulation implements not only continu- 
ous time equations but also discrete time equa- 
tions, thus the model has a mixed-data struc- 
ture. 

Introduction 

Read-time simulation of a realistic model is a 
cost effective way to design and test hardware. 
Model simulation is faster and safer than test- 
ing the actual system. Obviously, the com- 
puter system used to simulate a model is an 
important factor when considering simulation 
speed and simplicity pf implementation, but 
the simulation language is often equally im- 
portant. The ease with which a system can be 
modelled is strongly dependent on the simula- 
tion language. 

In this paper we present a realistic model us- 
ing the System 100. We try to show how con- 
veniently certain modelling problems can be 


handled using ADSIM. We discuss three major 
features of the aircraft model which are com- 
mon to many simulations. The first feature is 
the evaluation of aerodynamic functions with 
a large data base. This type of situation occurs 
in simulations where an analytic function does 
not exist and where function files dependent 
on several independent variables are given to 
describe a model. The second feature, the lon- 
gitudinal trim technique, is applicable to sys- 
tems where useful simulations must be run at 
an equilibrium condition. Finally, the case of 
mixed-data systems occurs whenever a digital 
and a continuous system are simulated. Al- 
though z-transform theory can be used to es- 
timate errors of a dynamic plant with a digital 
controller, it must be assumed that the plant 
is linear. If this is not the case, the plant must 
be approximated by a linear system. A thor- 
ough analysis of the error and system dynam- 
ics can be gained through simulation. Even if 
the plant is linear, oftentimes a knowledge of 
the inter-sample behavior is desired, and sim- 
ulation of the continuous plant and the digital 
controller is again required. 

In the next section a brief discussion of the 
System 100 is given. The following section 
gives an overview of the aircraft model; the ref- 
erence systems used, the orientation method, 
the external force model, and the control sys- 
tem model are each described. The last sec- 
tions are devoted to three specific features of 
the aircraft model, demonstrating how they 
are modelled using ADSIM. 
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The System 100 

The System 100 is an integrated simulation 
environment which consists of a general pur- 
pose computer and a digital simulation com- 
puter. The general purpose computer is one 
of the DEC VAX systems. The VAX front- 
end computer serves as a host for the Ap- 
plied Dynamics AD 100 and the ADSIM com- 
piler. Program preparation is done on the host 
computer. The AD 100 is the digital simu- 
lation compute engine and consists of up to 
seven parallel processors. The AD 100 is a 
totally synchronous, bus-oriented, multipro- 
cessor system capable of performing 20 mil- 
lion floating-point operations per second (20 
M Flops). 

The simulation language, ADSIM, is equation 
driven and block oriented. Many key elements 
of a typical simulation, such as integration 
techniques, function generation, and control 
system nonlinearity functions, are built into 
the language. A control executive consisting 
of two programs provides the basis for imple- 
menting a model. The two programs are IN- 
TEXEC which runs on the host computer and 
SIMEXEC which runs on the AD 100. The 
executive controls such parameters as simu- 
lation time, frame time, and integration step 
size. The user’s ADSIM program, consisting of 
blocks of code, is inserted into the SIMEXEC 
code at compile time. In this paper we men- 
tion two types of ADSIM blocks, the REGIOI 
block and the DYIAHIC block. REGIOI blocks 
may contain procedural code. The number of 
times and the order in which these blocks are 
executed are dependent on the specific name 
of the block. For instance, the code of REGIOI 
sync 2 is executed before the simulation run is 
begun while that of REGIOI sync4 is executed 
before each step of the simulation run. Eight 
optional REGIOI blocks are available. The 
DYIAHIC blocks contain the model dynamics. 
The code must be nonprocedural. The model 


differential equations reside here, written as a 
series of scalar and first order differential and 
difference equations. 

ADSIM offers an interactive environment 
which allows the program tote modified with- 
out recompiling. This environment is named 
IITERACT and allows instantaneous changes 
to simulation elements at run time which in- 
clude all numerical values of the program, in- 
tegration algorithms, integration step size, end 
time, sample time of a digital system and the 
speedup ratio with respect to real time. 

In this paragraph we note some ADSIM in- 
tegration terminology since we will need the 
concepts to describe the implementation of the 
model. These terms are referred to as run- 
specs of the ADSIM program([6]). The frame 
time is the amount of time needed for the 
computer to solve the differential model equa- 
tions. The step size, T, is the integration al- 
gorithm step used to calculate the system so- 
lution. The step time is the step size divided 
by the number of passes through the dynam- 
Tcal equations needed for the specific integra- 
tion algorithm. For instance the fourth order 
Runge-Kutta method, (RK-4), requires four 
passes through the differential equations and 
thus the step time is equal to ^ . In order for 
the simulation to run at real-time the actual 
time it takes to integrate the equations should 
be equal to the simulated time, thus the step 
time should equal the frame time. As an ex- 
ample, suppose the model uses the RK-4 inte- 
gration method, thus for real-time simulation 
the step size divided by four should be equal 
to the frame time. IF the step time is larger 
than the frame time the program is running 
at faster than real-time. 

Aircraft Model 

The aircraft model is representative of a busi- 
ness jet. The model can represent many other 
types of aircraft simply by changing some of 
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the simulation parameters. The jet aircraft is 
modelled as a rigid body, the center of mass 
chosen as the reference point so that the trans- 
lational and rotational aspects of the motion 
can be analyzed separately. The state vari- 
ables chosen to describe the translational mo- 
tion are: total velocity relative to the atmo- 
sphere, angle of attack, sideslip angle, distance 
north, distance east and altitude ([1]). The 
rotational motion is described using the con- 
ventional aircraft Euler angles and the com- 
ponents of the aircraft angular velocity along 
body axes. Figure 1 shows a block diagram of 
the overall six-degree-of-freedom flight equa- 
tions, Starting on the left are the state vari- 
ables and ending on the right are the time rates 
of change of the state variables. The time rates 
of change are integrated over each integration 
step to obtain the state variables used for the 
next step. 

A total of four reference systems is used in 
the model to conveniently describe the air- 
craft motion. The model orientation is best 
described using a Body reference system, a 
coordinate set fixed to the aircraft. A nonro- 
tating flat -earth atmospheric reference is de- 
fined to enable the translations and rotations 
of the aircraft to be related to an absolute ref- 
erence. The translational equations of motion 
are written using Flight Path axes since the 
total velocity, angle of attack, and sideslip an- 
gle are easily described in that frame. (The 
Flight Path frame is defined such that the to- 
tal aircraft velocity is along the X axis.) The 
Stability frame, an intermediate axis set be- 
tween the Body and the Flight Path frames, 
is used to describe the aerodynamic equations 
since force and moment data are usually ob- 
tained in that frame. 

Only three angles are needed to describe every 
possible orientation of the rigid body model. 
Of the choices of well known methods to 
describe rigid body orientation, this aircraft 
model implements the method of quaternions. 


The method of Euler angles was not used 
since a singularity occurs when the aircraft is 
pitched ±90 degrees. The method of direction 
cosines was not chosen since the method uses 
9 parameters to define the model orientation, 
requiring 6 constraint equations. The method 
of quaternions only requires 4 parameters to 
describe the system orientation, and thus only 
one constraint equation is needed. ADSIM 
offers a convenient quaternion model, which, 
given the initial Euler angles (defined for a 
3 — 2—1 rotation), and the aircraft’s present 
angular velocities, computes the present Euler 
angles and the direction cosines of the system. 
The external forces acting on the jet aircraft 
are modelled as a constant gravity force, air- 
craft thrust, and aerodynamic forces. Four- 
teen aerodynamic coefficients are used to de- 
scribe the aerodynamics of the aircraft. The 
coefficients are dependent on the aircraft an- 
gle of attack, sideslip angle, deflections of the 
control surfaces, (aileron, elevator, rudder and 
flaps), and mach number. 

The aircraft flight control system is described 
with a pitch, roll, and a yaw-damper control 
loop. Each control loop contains an actuator 
described by a proportional plus-rate feedback 
second-order system. In the case of the pitch 
and roll loops, a proportional plus rate feed- 
back is subtracted from the autopilot angle in- 
put to drive the controller. In the case of the 
yaw-damper loop, the yaw rate is sent through 
a high-pass filter and subtracted from the au- 
topilot rudder input to drive the system. 

Aerodynamic Function Interpolation 

Function evaluation is an important task for 
real-time simulation. In the aircraft model 
there is a large aerodynamic data base. Two 
variables, density and the reciprocal of the 
speed of sound, are dependent on only one 
variable, the aircraft altitude, and 14 aerody- 
namic coefficients, (such as the lift, side-force 
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and drag coefficients), are dependent on up 
to four independent variables. ADSIM pro- 
vides function generation capabilities for both 
equally and arbitrarily spaced functions on the 
AD 100. ADSIM allows multivariable func- 
tions with up to seven independent variables. 
Basically the way ADSIM evaluates a func- 
tion of some independent values is to first use 
a binary search technique to determine be- 
tween which given independent data points 
each input variable lies. Next interpolation 
techniques are used to evaluate the function 
at the given independent variable values. In 
order to perform function evaluation, indepen- 
dent variable and function tables must be de- 
clared. The aircraft simulation code contain- 
ing function declarations for the density, and 
drag and lift coefficients is shown below ([4]). 


(* Assign data files to 

breakpoint tables and functions *) 


FILES 
h_data 
a_data 
m_data 
de_data 
df _data 


*h_sixdof .bpt' , 

> a_s ixdof .bpt* , 
'm^sixdof .bpt' , 
'de.s ixdof .bpt * , 
’df_sixdof .bpt> 


FILES 

air_density 

lift_coeff 

drag_coeff 


* rho_s ixdof . fun 1 , 
’liftc.sixdof .fun 1 , 
'cd_s ixdof .fun' 


(* The dynamic block containing the 
continuous time equations *) 


(* Function generation breakpoint 
table definitions *) 


DYNAMIC Continuous 


INTERPOLATION_INTERVALS 

h_data( 129 OF 129) 
INTERPOLATION_INTERVALS 

a_data(33 OF 33) 
INTERPOLATION_INTERVALS 

m_data(65 OF 65) 
INTERPOLATION_INTERVALS 

de_data(3 OF 3) 
INTERPOL ATI ON_IHTERVALS 

df_data(3 OF 3) 


(* Function generation 

definitions *) 

INTERPOLATION.FUNCTIONS 

air_density(h_data) 
INTERPOLATION.FUNCTIONS 
lift_coeff (a_data,m_data > 
de_data , df _data) 
INTERPOLATION.FUNCTIONS 
drag_coef f (a_data, 

m_dat a , de.dat a , df _dat a) 


(* Evaluate air density *) 
rho = air.density(h) 


(* Aerodynamic coefficient 
evaluation *) 

cl = lift_coeff (a,m,de,df ) 

cd = drag_coeff (a,m,de,df ) 

END DYNAMIC continuous 


Longitudinal Trim 

The model contains an option to trim the air- 
craft in longitudinal flight. The aircraft is as- 
sumed to be in symmetric level flight. Longi- 
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tudinal trim is accomplished with the function, 
IEVTOI_RAP, supplied by ADSIM. For a given 
flight regime, (aircraft altitude and speed), 
IEWTON_RAP determines the angle of attack, 
thrust magnitude, and elevator deflection to 
cause the longitudinal accelerations a nd p itch 
moment to be zeroed. 

The function MEWTOI_RAP uses the Newton 
ftaphson method, a successive approximation 
process, to determine the trim values. This 
method has quadratic convergence. Although 
only a local method, the Newton Raphson al- 
gorithm works very well for the aircraft flight 
model, because the aircraft eq u at ions are well 
behaved and an inital guess for the trim values 
with some physical insight causes the method 
to converge within about four iterations. 

The function ¥EWTOH_RAP is called in the 
REGIOI sync2 block of ADSIM, thus the equi- 
librium values can be determined before the 
simulation run and the aircraft control is ini- 
tialized for a trimmed flight. 

Simulation of Mixed-Data System 

The si mulat ion of both a continuous and a 
digital system presents special modelling prob- 
lems, such as simulation of the computational 
delay times for the Analog to Digital (A/D) 
and Digital to Analog (D/A) conversions, and 
the simulation of sampling times. A special 
execute pair exists in ADSIM so that these 
modelling questions need not be re-solved for 
each model containing digital and continuous 
systems ([3]). The special execute pair named 
EXECUTE mixed_data controls the execution 
of the integration, the implementation of the 
delay times, and the order in which the dy- 
namic systems are solved. Note that the con- 
tinuous time equations should be solved on a 
frame-by-frame basis while the discrete time 
equations should be solved only every integer 
multiple of the frame time. This integer is 
dependent on the sampling rate of the digital 


system. 

The execute pair for mixed-data system simu- 
lation is used in ADSIM programs along with 
two dynamic blocks. Simulations of this type 
are easily partitioned into two or more sub- 
systems. One subsystem consists of algebraic 
and difference equations that represent the 
digital /discrete- time control laws. The other 
subsystem consists of the algebraic and differ- 
ential equations that describe the b ehavior of 
the continuous time system. When creating 
the control structure for simulation of mixed- 
data systems, the following factors were con- 
sidered; ' it: :: V 

• Numerical integ ration of disc ontinuities 
introduced by the digital controller resid- 
ing in the continuous time subsystem. 

• Selection of numerical integration 
method (s). 

• Adjustment of either the sample period or 
integration step ^ize. 

• Proper padding of simulation frames in 
order to maintain real-time. 

• Representation of computational delay in 
the digital subsystem. 

• Representation of errors introduced by 
A/D and D/A conversions. 

Let us discuss the selection of the integra- 
tion method. Single-step predictor type meth- 
ods which use the current and past deriva- 
tives to predict the new state values do so un- 
der the assumption of continuous derivatives. 
This assumption does not hold when simulat- 
ing mixed-data systems. Strong discontinuites 
are introduced by the control command as de- 
termined by the discrete time system. Self- 
starting methods, such as the classical, multi- 
step Runge-Kutta second, third and fourth 
order variety, do not rely on this assump- 
tion and thus these methods could be used. 
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However, the classical Runge-Kutta methods 
are not suited for real-time simulation since 
they require external inputs to the integra- 
tion method before they become available in 
time. Applied Dynamics developed Runge- 
Kutta Real-Time methods which have simi- 
lar accuracy and stability characteristics as 
the traditional versions but can be used for 
real-time simulations. Assume that we wish 
to solve the following vector differential equa- 
tion: 

X = f(x,u(t)) ( 1 ) 

Where x is the state vector and u is the input 
vector. The following difference equation de- 
fines the second order Runge-Kutta Real-Time 
(RKRT2) integration method: 

Xn+i = x n +Tf(x n+ $,u n+ i) (2) 

where 

T 

x n+ ±= x n + —f(x„,u„) (3) 

We note from these equations that when using 
the RKRT2 integration method for the sim- 
ulation of mixed data systems, we must as- 
sure the sample time is an integer multiple 
times the integration step size, T. The control 
structure of EXECUTE mixed_data will auto- 
matically adjust the sample time to insure this 
condition (only if the sample time is not an in- 
teger multiple of the frame time). 

We now consider the simulation of the digital 
controller’s computational delay. Actual mi- 
croprocessor hardware does not produce con- 
trol action to the continuous system instanta- 
neously. The most common method is to com- 
pute the control effort and output as soon as 
the results become available. However, when 
simulating delay time, the lag must be an inte- 
ger multiple of the continuous system integra- 
tion step time. If an estimate of the computa- 
tional delay for a particular piece of hardware 


is known, EXECUTE mixed_data will allow the 
effect of both fixed and time- varying computa- 
tional delays to be simulated on the continuous 
system model. The error introduced by the 
limited accuracy of the A/D and D/A conver- 
sions can be modelled using a nonlinear con- 
trol ADSIM function named Quantizer which 
generates a symmetric staircase function. 

The general format for ADSIM programs using 
the control structure of EXECUTE mixed_data 
follows: 

TITLE {Simulation Title} 

DYNAMIC discrete 

{ Difference and algebraic equations 

to represent digital controller.} 

END DYNAMIC discrete 

DYNAMIC continuous 

{ Differential and algebraic equations 
to represent continuous -time plant.} 

END DYNAMIC continuous 

EXECUTE 'mixed.data* 

We now describe the implementation of the 
digital controller in the aircraft model. In 
order to design the digital pitch control al- 
gorithm, knowledge of aircraft pitch behavior 
is needed. The longitudinal dynamics of air- 
craft are in general dominated by two phenom- 
ena; phugoid motion and short period pitch- 
ing motion. Phugoid motion is a very slow 
and lightly damped oscillatory motion in the 
vertical plane induced by mismatch between 
lift and drag of the airframe for a particu- 
lar flight condition. The short period pitch- 
ing motion occurs when the aircraft pitch an- 
gle is not aligned with the flight path. The 
aircraft dynamics will drive the pitch angle 
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back toward the flight path. The short period need not be represented in the ADSIM pro- 
pitching mode is much faster than the phugoid gram since the simulation control structure of 
mode. The continuous system pitch dynamics the Execute mixed_data pair supples this ef- 
were simulated and the frequency and damp- feet implicitly. The relevant ADSIM code is 
ing ratio of the two motions were estimated shown below, 
from the transients. It was found that the 

phugoid motion exhibits a natural frequency TITLE 'Six Degree-of -Freedom 
of about 0.012 Hertz with a damping ratio of Aircraft Simulation' 

0 . 1 , while the short period pitching mode fre- 
quency is about 0.45 Hertz with a damping COMMENT 
ratio of 0.35. The actuator is modelled by a 

proportional plus rate feedback second order This is an ADSIM program for simulating an 
system, its dynamics are characterized by a aircraft in six degrees of freedom, 
natural frequency of 5.0 Hertz and a damping The program runs on the Applied Dynamics 
ratio of 0.5. Using these estimates a transfer a D 100. Sixdof performs all the computations 
function representation was used to approxi- required for the six— degree— of —freedom 
mate the continuous system dynamics, where aircraft equations of motion, 
the input is the elevator deflection command 
and the output is the aircraft pitch angle. 

A proportional-integral-derivative (PID) type END COMMENT 
of digital controller is implemented in the air- 
craft simulation. The digital controller trans- 
fer function, D(z), can be defined using z- 

transform notation and has the following form, (* Equations of motion of the 


X a l + a 2 z 1 + a 3Z 2 

= l + M-' + M-* ' 


( 4 ) 


aircraft*) 
DYNAMIC Continuous 


The controller constants, a 2 , 03 , b x , b 2t 

were chosen using the transfer function ap- (* Flight control systems A 
proximation of the continuous system. These actuators: *) 

constants are easily refined at run-time using 

INTERACT. (* Elevator actuator *) 

As discussed in a previous section, the AD- 
SIM language supports _both ^continuous time dei = Deilim*SAT Cdeicom* 

derivatives and discrete time difference equa- Ideilim) 

tions. Therefore the D(z) transfer Tunction udec = Kde*(dei-des-Cdedot* 
is simply converted to ADSIM notation and dedot) 

placed into the discrete dynamic block. The dedot 9 = Udelira*SAT(udec* 

continuous dynamic block contains the equa-^ Iudelim) 

tions for the continuous model equations. The das' = dedot 

prime symbol denotes the differential opera- de = des+detrim 

tor, and the pound symbol denotes a shift 

in time, or a next-state variable. We note (* Roll control system *) 
that the transfer function for a zero order hold 


E 

I 
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daicom = Kphi*(phii-phi- 

Cp*p) 

dai « Dailim*SAT(daicom* 

Idailim) 

(* Aileron actuator *) 
udac=Kda* (dai-das-Cdadot*dadot) 


r 7 =lnv_lmlxzsqxlnv_lxxlz*(mz+ 
Ixzxlnv_lz*mx) 

END DYNAMIC Continuous 

DYNAMIC Discrete 

(* Pitch control system *) 

e_theta = -(thetai-theta 
-Cq*q) 

f_theta# = e_theta 
g_theta# = f _theta 

deicom = 

Bl*u_theta + B2*v_theta + 

Al*e_theta + A2*f_theta + 

A3*g_theta 

u_theta# = deicom 
v_theta# = u_theta 

END DYNAMIC Discrete 

EXECUTE , mixed_Data , 

Speed of Simulation 

The time needed for a single pass through the 
highly nonlinear dynamic equations is 138.1 
microseconds, corresponding to a frame rate of 
about 7000 frames per second. Since accurate 
real-time simulations of aircraft require about 


a 10 to 30 second frame rate([5]), it would be 
possible to run this model up to 700 times real 
time. 

Conclusions 

We have presented an ADSIM model for the 
real-time simulation of an aircraft. We have 
described the modelling of three specific fea- 
tures of the simulation which are common to 
many simulation models. These features are 
a large aerodynamic data base requiring func- 
tion interpolation, an option to trim the air- 
craft, and a digital controller. We have shown 
the model can be implemented and executed 
conveniently using ADSIM. 
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